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Non-modal amplification of disturbances in streamwise-constant channel flows of Oldroyd- 
B fluids is studied from an input-output point of view by analyzing the responses of the 
velocity components to spatio-temporal body forces. These inputs into the governing 
equations are assumed to be harmonic in the spanwise direction and stochastic in the 
wall-normal direction and in time. An explicit Reynolds number scaling of frequency 
responses from different forcing to different velocity components is developed, showing 
the same i?e-dependence as in Newtonian fluids. It is found that some of the frequency 
response components peak at non-zero temporal frequencies. This is in contrast to Newto- 
nian fluids, where peaks are always observed at zero frequency, suggesting that viscoelas- 
tic effects introduce additional timescales and promote development of flow patterns with 
smaller time constants than in Newtonian fluids. The temporal frequencies, correspond- 
ing to the peaks in the components of frequency response, decrease with an increase in 
viscosity ratio (ratio of solvent viscosity to total viscosity) and show maxima for non- 
zero elasticity number. Our analysis of the Reynolds- Orr equation demonstrates that the 
energy-exchange term involving the streamwise/wall-normal polymer stress component 
Txy and the wall-normal gradient of the streamwise velocity dyU becomes increasingly 
important relative to the Reynolds stress term as the elasticity number increases, and is 
thus the main driving force for amplification in flows with strong viscoelastic effects. 



1. Introduction 



Complex dynamical responses arise in numerous viscoelastic fluid flows ( Larson|[l992 



Shaqfeh 1996) and their study is important from both fundamental and technological 



standpoints. From the former standpoint, the inception and evolution of amplification 
of disturbances in various flows involving viscoelastic fluids is not well understood. Vis- 
coelastic effects not only modify features already present in Newtonian fluids, but also 
give rise to completely new behavioral patterns (Larson 2000 Groisman & Steinberg 



2000 Bertola et al. 2003). From the latter standpoint, the study of dynamics in flows 



involving polymeric fluids is of immense importance for polymer processing and rheom- 
etry (Bird et al. 1987 Larson 1999). Classical linear hydrodynamic stability analysis is 
found to give misleading results even for simple Couette and Poiseuille flows of Newto- 
nian fluids ( Schmid|[2007 ). This failure of the classical stability analysis is attributed to 



the non-normal nature of the generators in the linearized governing equations ( Trefethen 
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et al. 1993 'Grossmann""2000 Schmid & Henningson 2001). Linear dynamical systems 



with non- normal generators can have solutions that grow substantially at short times, 
even though they decay at long times ( Gustavsson||1991 Butler fc Farrell||1992 Reddy 



fc Henningson 1993j). Furthermore, the non-normal nature of the underlying equations 
can lead to significant amplification of ambient disturbances ( Farrell fc loannou 1993| 



Bamieh fc Dahleh]|2001 Jovanovic fc Bamieh||2005 ) and substantial decrease of stability 
margms dTrefethen et a/.|1993| [Trefethen fc Embree|2005D . On some occasions, the tran- 
sient growth and amplification, which are overlooked in standard linear stability analysis, 
could put the system in a regime where nonlinear interactions are no longer negligible. 
These phenomena are also expected to be important in Couette and Poiseuille flows of 
viscoelastic fluids. In this paper, we investigate amplification of disturbances in channel 
flows of Oldroyd-B fluids by performing a frequency-response analysis. 

Novel ways of describing fluid stability that allow quantitative description of short-time 
behavior and disturbance amplification, referred to as non-modal stability analysis, have 
emerged in the last decade ( Schmid|2007[ ). One approach is to study the responses of the 
linearized Navier-Stokes equations (LNSE) to external disturbances (Farrell fc loannou 



1993 Bamieh fc Dahleh 200T] ). Jovanovic fc Bamieh ( 2005) have used this approach to 



study the effects of external disturbances, in the form of body forces, on channel flows 
of Newtonian fluids. Explicit Reynolds number dependence of the components of the 
frequency response was derived. Based on this finding, it was concluded that at higher 
Reynolds numbers, wall-normal and spanwise disturbances have the strongest influence 
on the flow field and the impact of these forces is largest on the streamwise velocity. It 
was also found that the frequency responses peak at different locations in the ikx^kz)- 
plane, where kx and kz are the streamwise and spanwise wavenumbers, indicating the 
possibility of distinct amplification mechanisms. We note that - even in high- Reynolds- 
number regimes - it is valid to examine the linearized equations to determine the fate of 
small-amplitude perturbations to the underlying base flow. 

Recently, the present authors have extended the work of Jovanovic fc Bamieh" ("2005") to 
viscoelastic fluids (Hoda, Jovanovic fc Kumar 2008). Prior studies on transient growth 
phenomena in viscoelastic fluids were reviewed there, and hence will not be reviewed 
here for brevity. The aggregate effect of stochastic disturbances in all the three spatial 
directions to all the three velocity components, referred to as the ensemble-average en- 
ergy density, was investigated. It was found that the energy density increases with an 
increase in elasticity number and a decrease in viscosity ratio (ratio of solvent viscosity to 
total viscosity). In most of the cases, streamwise-constant or nearly streamwise-constant 
perturbations are most amplified and the location of maximum energy density shifts 
to higher spanwise wavenumbers with an increase in elasticity number and a decrease 
in viscosity ratio. However, prior work on Newtonian fluids by [Jovanovic fc Bamieh 



( 2005 ) suggests that a plethora of additional insight can be uncovered by analyzing the 
componentwise spatio-temporal frequency responses. The componentwise responses give 
information about the relative importance of the three disturbances on the three ve- 
locity components. By analyzing these frequency responses, the disturbance frequency 
corresponding to maximum amplification can also be obtained. 

In this paper, we examine the componentwise frequency responses in streamwise- 
constant channel flows of Oldroyd-B fluids. This study supplements a previous study 
by the authors (Hoda et al. 2008"), where the aggregate effect of disturbances was ex- 
amined, and helps in understanding the relative importance of the disturbances on the 
different velocity components. Because the previous study focused on aggregate effects, as 
parameterized by the energy density, it leaves open the question of exactly which velocity 
components are most amplified and which forcing components are responsible for this 
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v, = 2U„[l-(y/Lf] 



Couette Flow: 
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Figure 1 . Schematic of the channel- flow geometry. 



amplification. Furthermore, since the energy density is a time-integrated quantity, it does 
not yield information about most amplified temporal frequencies. The present work ad- 
dresses these issues, provides some explicit scaling relationships, and further investigates 
physical mechanisms. 

Streamwise-constant three-dimensional perturbations are considered in this work as 
they are most amplified by the linearized dynamics. We derive an explicit Reynolds 
number scaling for the components of the frequency response. As in Newtonian fluids, 
at higher Reynolds numbers the forces in the wall-normal and spanwise directions have 
the strongest influence on the flow field and the impact of these forces is largest on 
the streamwise velocity. In some of the cases, the frequency response components peak 
at non-zero temporal frequencies. This is distinct from Newtonian fluids, where peaks 
are always observed at zero frequency, suggesting that elasticity introduces additional 
timescales and promotes development of flow patterns with smaller time constants than 
in Newtonian fluids. We also find that the temporal frequencies, corresponding to the 
peaks in the components of frequency response, decrease with an increase in viscosity 
ratio and show maxima with respect to the elasticity number. One of the most important 
conclusions of this paper is the observation that elasticity can lead to considerable energy 
amplification even when inertial effects are weak; this energy amplification may then serve 
as a route through which channel flows of Oldroyd-B fluids transition to turbulence at 
low Reynolds numbers ( Larson||2QQQ Groisman fc St einberg||2QQQ| . 

Our presentation is organized as follows: in §[2j a model for streamwise-constant chan- 
nel flows of Oldroyd-B fluids with external forcing is presented. In § [3j a brief summary 
of the notion of the spatio-temporal frequency response is provided. In § [4j an explicit 
scaling of the frequency response components with the Reynolds number is given. In 
§ |5] and § |6j the effects of elasticity number and viscosity ratio on power spectral and 
steady-state energy densities are studied. The important findings are summarized in § [7| 
and the detailed mathematical derivations are relegated to the Appendix. 



2. The cross-sectional 2D/3C model 

A schematic of the channel-flow geometry is shown in figure [T] The height of the 
channel is 2L and the channel extends infinitely in the x- and z- directions. For Couette 
flow, the top plate moves in the positive x-direction and the bottom plate moves in 
the negative x-direction, each with uniform velocity Uq- For Poiseuille flow, 2Uo is the 
centerline velocity. We analyze the dynamical properties of the LNSE for an Oldroyd-B 
fluid with spatially distributed and temporally varying body-force fields. The parameters 
characterizing Oldroyd-B fluids are: (a) the viscosity ratio, (5 — rjs/ (Vs-^Vp)^ where r]s and 
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r]p are the solvent and polymer viscosities, respectively; (b) the Reynolds number, Re = 
UoL/v = pUoL/ {rjs -^Tjp)^ which represents the ratio of inertial to viscous forces, where p 
denotes fluid density; (c) the Weissenberg number. We = XUo/L^ which characterizes the 
importance of the fluid relaxation time. A, with respect to the characteristic flow time, 
L/Uq. Another important parameter is the elasticity number, p = We /Re = Az^/L^, 
which quantifies the ratio between the fluid relaxation time. A, and the vorticity diffusion 
time, L'^/u. 

Consider the dimensionless linearized momentum, continuity, and constitutive equa- 
tions 

-v-Vv - v-Vv + ^f-Vp 
_ Re\ 
V -V, 

^ - (Vv + (Vv)^) - V- Vr 

r 

We' 

where v = [u v w ]^ is the velocity fluctuation vector, p is the pressure fluctuation, 
and r is the polymer stress fluctuation. The overbar denotes the base flow given by 

zvve [U [V}^'^ 



dtv = 
= 

dtT = 



We 

(r- Vv)^ + (r - Vv)^ 



/^V^v + (l-/3)V-r) + d. 



V - Vr + T - Vv + T - Vv 



(2.1) 



V = [ U{y) Y 



T = 



2We{U\y)Y 




U\y) 






with U{y) = {y, Couette flow; 1 — y^^ Poiseuille flow}, and U'{y) = dU{y)/dy. A 
spatio-temporal body force is represented by d = [ ^2 ds ] , where di,d2, and ds 
are the body force fluctuations in the streamwise (x), wall- normal (y), and spanwise (z) 
directions, respectively. These body forces can be either deterministic or stochastic, and 
they serve as inputs into the system of equations that governs evolution of velocity and 



polymer stress fluctuations (Jovanovic & Bamieh 2005^ . Our objective is to investigate 
their effect on the components of the velocity field ^r = [ u v w ]^ . 

In this paper, we confine our analysis to streamwise-constant perturbations. In this spe- 
cial case, the model for flow perturbations is usually referred to as the two-dimensional, 
three-component (2D/3C) model (Reynolds & Kassinos 1995[ ). (2D indicates that the 
dynamics evolve in the cross-sectional {y, z)-plane, and 3C indicates that velocity com- 
ponents in all three spatial directions are considered.) The motivation for a thorough 
analysis of this particular model is twofold: (a) a recent study by the authors suggests that 
streamwise-constant (and nearly streamwise-constant) perturbations in channel flows of 
Oldroyd-B fluids create the largest contribution to the ensemble- average energy den- 
sity ( Hoda et a/.||2QQ8 ); and (b) for the 2D/3C model, an explicit i?e-dependence for the 
components of the frequency response can be obtained, which clarifies the effectiveness 
(energy content) of forcing (velocity) components. 



The evolution model of the forced linearized system (2.1) with streamwise-constant 
perturbations {kx = 0) is obtained by a standard conversion ( [Schmid fc Henningso n 
|2QQ1| ) to the wall-normal velocity /vorticity {v^cOy) formulation: 



dt^l^{y,kz,t) = 
v{y,k^,t) = 



[A{k,)^l,{k,,t)]{y) ^ 
[C{k,)tl^,{k,,t)]{y). 



[B{k,)d{k,,t)]{y), 



(2.2) 



Here, kz is the spanwise wavenumber, i/? = 



1 T 



^ il)i = [v ujy] , and the com- 
ponents of polymer stress are given by i/?2 = [^xx ^yy ^zz ^xy ^xz ^yz]^ • The derivation 



of the evolution equation for v requires elimination of the pressure from Eqs. (2.1). This 
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is achieved by applying the divergence operator to the momentum equation and by com- 
bining the resulting equation with continuity. On the other hand, the equation for cOy is 
derived by applying the curl operator to the momentum equation. 



The operator A in Eq. (2.2) is referred to as the dynamical generator of the linearized 
dynamics and it characterizes internal properties of the LNSE (e.g., modal stability). 
The definition of this operator for full three-dimensional fluctuations is provided in |Hoda| 
et al. ( 2008 ) ; the deflnition of components of this operator suitable for frequency response 



analysis of the 2D/3C model is given in Eqs. (2.3) and (2.4) below. We also note that 
operator B can be partitioned as B = [ B^ ] , where B describes how forcing 
enters into the Orr-Sommerfeld and Squire equations of viscoelastic channel flows, and 



O is a 6 X 3 matrix of null operators. On the other hand, operator C in Eq. (2.2 ) contains 
information about a kinematic relationship between '0^ and v. These two operators are 
given by: 





Bi 



B2 

























{Bi = B2 = -klA-\ B3 = -ifc.A-i^j,}, {C„ = -(i/fc,), C„ = I, C» = {i/k,)dy}, 

where i = a/— T, I is the identity operator, A = dyy — /c^ is a Laplacian wi th D irichlet 
boundary conditions, and denotes the inverse of the Laplacian. System (2.2) is sub- 



ject to the following boundary conditions {'^(±1, kz^t) = dyv{±.l^ kz,t) = cjy(±l, kz^t) = 
0}, which come from the no-slip and no-penetration requirements. We note that no 
boundary conditions on the polymer stresses are needed ( Hoda et a/. [2008 ). 

A coordinate transformation = Ti/? with {^1 = 02 = l^yy ^yz ^zz]^ ? ^3 = 
04 = i^xy Txz ] , ^5 = 'Txx} cau be uscd to bring system (2.2) into the following form: 



u 

V 

w 



Re 



Re 





4>i + 


1 




jj^Re 


02 + 


- F33 


h + 


1 




jiiRe 


04 + 


1 




jj^Re 


4>5 - 

















;iven 


by: 




F33 



Re 
1 



Fl2 02 



B2d2 



Bsds, 



jiRe 
F31 01 



F21 01, 
1 



/5 



'41 ^1 



fiReY^i ( 



Re 

F42 02 



F34 04 + Bidi 



1 



M3 ^3, 



jiRe 

F53 03 + F54 04, 



(2.3a) 
(2.36) 
(2.3c) 
{2.U) 
(2.3e) 

(2.3/) 



Fn = A-^A^ F33 = A, F31 = -ikzU'{y), 

F12 = A-l [ -kldy -ikz {dyy + k^) k^dy ] , F34 = [ 1^.^^ ' k^ ] 

F2I = [ 2dy {i/kz){dyy + k^) -2dy F43 = ' {l / k^Yl,, 



(2.4) 



U\y)dy - U"{y) 

(\/kz)U'{y)dyy 

F51 = W{y)U"{y), F53 = - 



M2 



U'{y) 
U\y) 
(2\/kz)U'{y)dy, F54 = [ 2U'{y) ] . 
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Here, A" 



yyyy 



- with both Dirichlet and Neumann boundary conditions. 
The system of equations (2.3) is in a form suitable for the analysis performed in § |4] 



yy 



where an explicit characterization of the Reynolds-number dependence for the compo- 



nents of the frequency response of system (2.2) is provided. It is noteworthy that for 
the 2D/3C model there is no coupling from ^5 = Txx to the equations for the other 
flow-field components in (2.3); in particular, this demonstrates that evolution of Txx at 



kx = does not infiuence evolution of ii, and w. We also note a one-way coupling from 
Eqs. (2.3a) and (2.3b) to Eqs. (2.3c) and (2.3d); this indicates that the dynamical prop- 

T 



erties of ^3 = 
but not vice- versa. 



and 04 = [ Txy Txz ] are infiuenced by ^1 = v and 02 = [ Tyy Tyz r^z ] 



3. Frequency responses for streamwise-constant perturbations 

Frequency response represents a cornerstone of input-output analysis of linear dynam- 



ical systems (Zhou et al 1996). The utility of input-output analysis in understanding 



early stages of transition in wall-bounded shear fiows of Newtonian fiuids is by now well 



documented; we refer the reader to a recent review article by Schmid (2007) for more 



information. It turns out that the input-output approach also reveals important facets 
of transitional dynamics in channel fiows of Oldroyd-B fiuids ( Hoda et a/.||2008 ). 

To provide a self-contained treatment, we next present a brief summary of the notion 
of the spatio-temporal frequency response of the streamwise-constant LNSE with forcing; 
we invite the reader to see Jovanovic & Bamieh ( 2005| for additional details. The spatio- 
temporal frequency response of system (|2.2|) is given by 



H(fc„c^) = C{k,){iiol-A{k,))-'B{k,), 
where uj denotes the temporal frequency. The frequency response is obtained directly from 



the Fourier symbols of the operators in Eq. (2.2), and for any pair (kz^oo) it represents 
an operator (in y) that maps the forcing field into the velocity field. 

The frequency response of a system with a stable generator A describes the steady- 
state response to harmonic input signals across temporal and spatial frequencies. Since 
H is an operator valued function of two independent variables kz and cj, there are a 
variety of ways to visualize its properties. In this paper, we study the Hilbert-Schmidt 
norm of H 

U{kz,uj) = trace (H(A:^,cj)H*(A:^,cc;)), (3.1) 

where H* represents the adjoint of operator H. For any pair (/c^, cj), the Hilbert-Schmidt 
norm quantifies the power spectral density of the velocity field in the LNSE subject to 
harmonic (in z) white, unit variance, temporally stationary, stochastic (in y and t) body 
forcing. Furthermore, the temporal-average of the power spectral density of H yields the 
so-called H2 norm of system (2.2) ( |Zhou "erZ][T996l ) 



)duj. 



The frequency responses of viscoelastic channel fiows (as a function of kx and kz) are 
quantified in Hoda et al. (2008) in terms of the H2 norm. We note that at any /c^, the 
H2 norm determines the energy (variance) amplification of harmonic (in z) stochastic 
(in y and t) disturbances ( Farrell fc Ioannou| 1993 Bamieh fc Dahleh||2001[ |Jovanovi6 
|fc Bamie"h||2005| ). This quantity is also known as the ensemble- average energy density of 
the statistical steady-state ^ Farrell fc Ioannou|1993 ), and it is hereafter referred to as the 
(steady-state) energy density (or energy amplification). 
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We finally note that the frequency response of system (2.2), v 
3x3 block-decomposition: 



Hd, has the following 



u 

V 

w 



^ui(kz,uj] Re, 13,11) 
^yi(kz,u]Re,(3,ii) 



^u2(kz,uj',Re,l3,ii) 
Hy2{kz,uj;Re,(3,iJ.) 



Hn3(^2,^;^e,/3,/i) 
lly3{kz,u;;Re,p,fi) 













. ^3 



(3.2) 

which is suitable for uncovering the effectiveness (energy content) of forcing (velocity) 
components. In this representation, H^j denotes the frequency response operator from 
dj to r, with {j = 1,2,3; r = u,v,w}. Our notation suggests that — in addition to the 
spanwise wavenumber kz and the temporal frequency uj — each component of H also 
depends on the Reynolds number Re, the viscosity ratio P, and the elasticity number fi. 



4. Dependence of frequency responses on the Reynolds number 

In this section, we study how the power spectral densities and the steady-state energy 
densities scale with Re for each of the components of the frequency response (3.2). Fur- 
thermore, the square- additive property of these two quantities is used to determine the 
aggregate effect of forces in all three spatial directions d = [ di ^2 c^s ]^ on all three 
velocity components v = [ u v w ]^ . We analytically establish that the frequency re- 
sponses from both wall- normal and spanwise forces to streamwise velocity scale as Re^, 
while the frequency responses of all other components in Eq. (3.2) scale at most as Re. 
This extends the Newtonian- fluid results ( Jovanovi6||2QQ4 Jovanovic fc Bamieh||2QQ5 ) to 
channel flows of Oldroyd-B fluids. 

Application of the temporal Fourier transform to Eq. (2.3) facilitates elimination of 
polymer stresses from the 2D/3C model (see Appendix A| for details). This leads to 
an equivalent representation of system ( |2.3[ ) in terms of its block diagram, which is 
shown in figure [2] with Q = uoRe. From this block diagram, it follows that operator 
'tl{kz,uj; Re, P, /i) in Eq. (|3.2[) can be expressed as 



u 




Re llui{kz,^',P,ii) 


Re'^]iu2{kz,n',P,li) 


Re'^]iu3{kz,n',p,/i) " 




di 


V 







Re 'H.y2{kz,^', (3, fi) 


Re liys{kz,^;P,ii) 




d2 


w 







Rei:iy,2{kz,Vt](3,ii) 


Rellyjs{kz,ft;P,ii) 




ds 



where the Reynolds-number-independent operators H are given by: 
H„i(fc„!^;/3,//) = C„(im - S)-iBi, 



(4.1) 



H'fxj (^kz 



Urj{kz,n;p,fi) = Cr{ini - 
Operators L, S, and Gp are defined by: 
1 + _i . 9 ^ 1 



L)-iB„ {r 



1 + ifift 



-ikzU'iv), Cp2 



1 



f i/iQ 
ikz{l 



P){U'{y)A 



i = 2,3, 
2,3}. 

A* 

2U"{y)dy) . 



In the limit p ^ 1, these operators simplify to the familiar Orr-Sommerfeld (A~^A^), 
Squire (A), and coupling {—\kzU'{y)) operators in the streamwise-constant LNSE of 
Newtonian fluids with Re = 1. We note that even though viscoelastic effects modify 
some of the operators in figure [2j there is a striking similarity between block-diagram 
representations of the 2D/3C models of non-Newtonian and Newtonian fluids (jJovanovic 



fc Bamieh||2QQ5 ). In particular, figure [2] shows that the frequency responses from d2 and 



Bi 
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C!2 


B2 


>H 










d3 


B3 










Figure 2. Block diagram representation of the streamwise constant LNS system (2.2) 



ds to u scale as i?e^, whereas the responses from all other forcing components to other ve- 
locity components scale linearly with Re. It should be noted that the frequency responses 
of Newtonian fluids =0 show same scaling with Re (| Jovanovic || 2 004 Jovanovic & 

Bamieh||2005[ ). The coupling operator, Cp, is crucial for the i?e^-scaling. In Newtonian 
fluids Cp corresponds to the vortex tilting term, Cpi; in viscoelastic fluids Cp also con- 
tains an additional term, Cp2, that captures the coupling from the wall- normal velocity 
to the wall-normal vorticity due to the work done by the polymer stresses on the flow. 
In the absence of the coupling operator, all the components of H(cc;, /c^; i^e, /i) scale 
at most linearly with Re. Figure [2] also suggests that for the 2D/3C model, streamwise 
forcing does not influence the wall-normal and spanwise velocities, which is in agree- 
ment with Newtonian-fluid results ( [Jovanovic fc Bamieh|2005 ). As noted in Appendix [Aj 
Gp2 has its origin in the term involving polymer stress fluctuations and gradients in the 
base velocity proflle. This would produce polymer stretching, and could be interpreted 
as giving rise to an effective lift-up mechanism (Landahl||1975| ). 

The i?e-scaling for the power spectral densities of operators H^j(/c2, cj; i?e, /i) in 
streamwise-constant Poiseuille and Couette flows of Oldroyd-B fluids follows directly 



from Eqs. (3.1) and (4.1) and linearity of the trace operator: 



Uui{kz,uj;Re,(3,iJ.) 
Uyi{kz,uj]Re,p,/j.) 
U^i{kz,u;;Re,p,iii) 

U^i{k,,n;P,fi)Re^ 





Uu2{kz,uj;Re,(3,iJ.) 
Uy2{kz,uj]Re,p,/j.) 

U^2{kz,n;P,^)Re^ 
Uy2{kz,^;f^,l^)Re'^ 
Uu,2{kz,ft;p,/j.)Re'^ 



Ilu3{kz,uj;Re,f3,jj.) 
Uys{kz,uj;Re,p,/j.) 
U^3{kz,uj;Re,p,/i) 

U^s{kzM;P.f^)Re^ 
Uys{kz,Q;P,lii)Re'^ 
Uu,3{kz,n;p,/j.)Re'^ 



(4.2) 



where ILrj are the power spectral densities of the Reynolds-number-independent opera- 
tors 'H.rj{kz^ ^; /i), with Q = uoRe. Furthermore, the power spectral density of operator 
H(/c2, cj; i?e, /3, /i), v = Hd, is given by 

n(fe^,cj;i?e,A/i) = Iia{kz,9^]P,^)Re^ + n5(fc^, 1^; A /i)i?e^ 



where 11^ = II^i 



n 



v2 



n 



v3 



n 



w2 



n^3 and 



n 



n3 • 



Several important observations can be made about Eq. (4.2) without doing any de- 
tailed calculations. First, the power spectral densities of operators H^2 and H^3 scale as 
Re^] in all other cases they scale at most as Re^ . This illustrates the dominance of the 
streamwise velocity perturbations and the forces in the wall-normal and spanwise direc- 
tions in high-Reynolds-number channel flows of streamwise constant Oldroyd-B fluids. 
Second, apart from Iiu2{kz^ Vt] /3, /i) and Iiu?,{kz^ O.] /i), the other power spectral densi- 
ties in Eq. (4.2 ) do not depend on the base velocity and stresses. These two power spectral 



Frequency responses in streamwise- constant channel flows of Oldroyd-B fluids 9 



densities depend on the coupling operator, Cp, and thus their values differ in Poiseuille 
and Couette flows. Third, power spectral densities ILrj{kz^^] ja) do not depend on 
the Reynolds number. Thus, Re only affects the magnitudes of lirjikz^oo^Re^ jS^ ji)^ and 
regions of temporal frequencies uo where these power spectral densities peak. As Re in- 
creases, these C(;-regions shrink as 1/Re. Therefore, for high- Reynolds- number channel 
flows of Oldroyd-B fluids, the influence of small temporal frequencies dominates the evo- 
lution of the velocity perturbations, suggesting preeminence of the effects in fluids with 
relatively large time constants. It is noteworthy that elasticity shifts temporal frequencies 
where fl^j peak to higher values, which is discussed in detail in §[5] For additional details 
concerning these points, we refer the reader to Hodaj ([2008j) . 

We next exploit the above results to establish the Reynolds- number dependence of 



steady-state energy densities for different components of frequency response operator (3.2 ) 
For example, £^^2(^2; ^e, /i) is determined by 



Eu2{kz]Re,l3,ii) 



1 

— j Uu2{uJ,kz;Re,p,/j.)du; 

/CO 
-CXD 



27T 

Re' 



-00 

_ Re^ 
=: Re^gu2{kz',P,lj). 

A similar procedure can be used to determine the steady-state energy densities of all 
other components of operator H in Eq. ( |3.2[ ), which yields: 

Euiikz] Re, 11) Eu2{kz;Re, fi) E^sikz; Re, /3, /i) 
Eyi{kz; Re, f5, fi) Ey2{kz]Re, p, fi) Ey^{kz] Re, p, /i) 
Ewi{kz', Re, p, 11) Ey,2{kz] Re, p, /i) E^s{kz', Re, p, /i) 

fui {kz] P, fi)Re gu2{kz; P, /^)Re^ 9u?>{kz] P, fi)Re^ 
fy2{kz;P,fi)Re fy^{kz', P, ij)Re 

fw2{kz',P,/J^)Re U3{kz',P,/J^)Re 

where frj and grj are functions independent of Re. Furthermore, the steady-state energy 
density of operator H{kz,uj; Re, P, /j.), v = Hd, is given by 



E{kz;Re,P,fi) = f{kz;P,fi)Re + g{kz; P, fi)Re^ , 



(4.3) 



where / = /ni + fv2 + fvs + fw2 + fws and g = gu2 ^ QuS- 

We conclude that energy amplification from both spanwise and wall-normal forcing 
to streamwise velocity is 0{Re^), while energy amplification for all other components is 
0{Re). 



5. Parametric study of power spectral densities 

In §|4j we derived an explicit dependence for each component of the frequency response 
operator ( |3.2[ ) on the Reynolds number. Here, we investigate the effect of p and /i on 
the (O, /c^j-parameterized plots of power spectral densities fl^j, {r = u, v, w, j = 1, 2, 3}, 
by setting Re = 1 in Eq. (4.2). In all the plots presented in this section, 100 x 90 log- 
arithmically spaced grid points are used in the (1], /c^) -plane. The temporal frequency 
and spanwise wavenumber are varied between 0.01 and 25.11 (ft) and 0.1 and 15.84 (kz), 
respectively. The Reynolds- number- independent power spectral densities fLrj{kz, P, j^) 



10 



N. Hoda, M. R. Jovanovic, and S. Kumar 



in Eq. (4.2) can either be numerically determined from a finite-dimensional approxima- 
tion of the underlying operators or they can be computed using the method developed 
by Jovanovic & Bamieh (2006). For numerical approximation, we use a Chebyshev collo- 
cation technique (Weideman & Reddy 2000); between 30 and 50 collocation points were 
found to be sufficient to obtain accurate results. In Couette flow, the power spectral 
densities can be computed more efficiently using the method developed by [Jovanovic fc| 



Bamieh 



(|2006l). This can be accomplished by expressing each component of the frequency 
response op erator in a d ifferent form, known as the two-point boundary value state-space 
realization (|Hoda|2008). In Couette flow, we used both methods for evaluating 11^^; the 



results agreed with each other, suggesting accuracy of our computations. 

Figures [3] and |4j respectivel y, sho w the (1], /c^j-dependence of the i?e-independent power 
spectral densities ILrj in Eq. (4.2), for P = 0.1 and /i = 10. As noted in § [ij only n^2 
and n^3 depend on the base velocity and polymer stresses. In view of this, the results 
in figure |4] are computed in both Couette and Poiseuille flows. Since, at kx = 0, di does 
not affect v and we do not plot fl^i and fl^i in figure [s] Also, since n^2 = H^a, we 
only plot Ilys{kz^^; /J^). We now discuss some important observations concerning the 
results presented in these figures. 

It is clearly seen that several frequency-response components peak at non-zero Q val- 
ues. This is in contrast to Newtonian fluids, where all power spectral densities attain 
their respective maxima at = ( Jovanovic fc Bamieh 2006) . Also, since the peaks 
for different components of the frequency response are observed at different locations 
in the (1^, /c^) -plane, these plots suggest distinct amplification mechanisms. It is worth 
mentioning that the locations of the peaks shift depending upon the /3 and /i values. Our 
results indicate that viscoelastic effects introduce additional timescales which promote 
development of spatio-temporal flow patterns with smaller time constants compared to 
Newtonian fluids. 

At low Re values (^1), depending upon /i, either input-output amplification from di to 
u attains the largest value or the amplification from (6/2, ds) to u attains the largest value. 
At small values of /i, fl^i has the largest magnitude; this suggests that at small Reynolds 
numbers and small elasticity numbers, the streamwise forcing has the strongest influence 
(on the velocity) and the most powerful impact of this forcing is on the streamwise 
velocity component. At higher values of /i, n^2 and ILus achieve the largest magnitudes; 
this suggests that at higher elasticity numbers, the spanwise and wall- normal forces have 
the strongest influence (on the velocity) and that the streamwise velocity component is 
most energetic. 

As discussed above, the power spectral density peaks (for different components of the 
frequency response operator) are observed at different locations in the (1], /c2)-plane. We 
next analyze the effects of parameters fi and P on the magnitude of these peaks and the 
locations of the respective maxima. Figures |5] and [6j respectively, illustrate the variation 
with /i (for p — {0, 1, 0.5, 0.9}) in the 1^- value corresponding to the maxima of functions 
ILrj in Eq. (4.2); we denote this value by l^max- Below, we discuss the key features of 
these results. 

For fi greater than a certain threshold value, l^max exhibits a maximum in /i. Also, the 
value of fi corresponding to the maximum in f^max decreases with a decrease in p. The 
above results suggest that for a particular range of elasticity numbers, viscoelastic effects 
in Oldroyd-B fluids promote amplification of flow structures with smaller time constants 
than in Newtonian fluids. For small elasticity numbers, figures [5] and [6] demonstrate that 
all power spectral densities achieve their peak values at zero temporal frequency. This 
is consistent with the behavior of Newtonian fluids (Jovanovic & Bamieh 2006) as an 
Oldroyd-B fluid is equivalent to a Newtonian fluid in the limit /i ^ 0. The plots in 
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Figure 3. Plots of the base-flow-independent power spectral densities Urj in Eq. (4.2) for 

(3 = 0.1, /i = 10. 



figure [5] also suggest that for large elasticity numbers, l^max monotonically decreases 
(with a slow rate of decay) as ja increases. 

It is also seen that l^max increases with a decrease in suggesting the importance of 
effects with shorter time constants in viscoelastic fluids. In the limit /3 ^ 1, l^max = 



(results not shown), which is in agreement with the behavior of Newtonian fluids (Jo- 



|vanovic fc Bamieh 2006); an Oldroyd-B fluid is equivalent to a Newtonian fluid in the 
limit p 1. 

For the base-flow-independent power spectral densities in Eq. ( |4.2[ ), a very good ana- 
lytical estimate for l^max can be determined by projecting the operators in H^j on the 

(for n^j, {r = v^w; j = 2,3}) and A (for II^i). (We refer 



first eigenfunctions of A ^A^ 



the reader to Appendix B of|Jovanovic fc Bamieh (2005) for spectral analysis of these 



two operators in the 2D/3C model.) Using this approach, we determine the following 
expression for l^max 



v//i|Ai(fe,)|(l - p){fi\Xi{h)\{l +^+2) - 1 
0, 



2 

1-/3 



1 



^ ^ |Ai(fe,)|(l + /?)' 
otherwise, 



where Ai(A:;^) is the principal eigenvalue of the underlying operator (A for II^^i; A ^A^ 
for n^j, {r = v^w; j = 2,3}). From this expression, it follows that l^max increases with 
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Figure 4. Plots of the base-flow-dependent power spectral densities Ilu2 and n^^s in Eq. (4.2) 
for Couette (first row) and Poiseuille (second row) flows with p = 0.1, /x = 10. 



a decrease in (3 and that it exhibits a maximum in /i, explaining the trends observed in 
figure|5] Furthermore, for /a ^ 1, ^max approximately scales as l/yT^, which justifies our 
earlier claim regarding the slow rate of decay of l^max with /i for large elasticity numbers. 

The analytical expressions for l^max are much more difficult to obtain for the base-flow- 
dependent power spectral densities n^2 and n^3. In spite of this, essential trends can be 
ascertained from figure [7| which shows the plots of function Il}j{kz,^; /n) in Couette 
flow with p = 0.1 and fi = {0.1, 2, 10, 100}. (This function quantifies the power spectral 
density of the frequency response operator that maps (^2, ds) to streamwise velocity u at 
Re = 1.) For /i = 0.1, the frequency response achieves a global maximum at l^max = 0, 
but the broad spectrum in Q around kz ~ 0(1) indicates that the large values are 
maintained up until ft ^ 5. On the other hand, at /i = 2, the global peak is located in 
the narrow region around l^max ^ 2. With a further increase in elasticity number, two 
competing peaks at zero and 0(1) temporal frequencies appear; finally, for large values 
of /i the spectrum peak shifts to the narrow region around zero temporal frequency, with 

^max — 0- 

We have also studied the effects of ja and p on the maximum value of power spectral 
densities. We briefly highlight several important points. First, the maximum values of all 
power spectral densities decrease with an increase in P and a decrease in /i. This suggests 
that amplification becomes weaker as one approaches the Newtonian fluid limit. A simi- 
lar dependence of the maximum value of the growth function on Deborah (Weissenberg) 
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Figure 5. Variation in Qmax with /i for /3 = {0.1, 0.5, 0.9}. 



number was reported by Sureshkumar et al. (1999) in their study of two-dimensional 



time-dependent simulations of creeping plane Couette flow of Oldroyd-B fluids. Second, 
the peak values of the base-flow-dependent power spectral densities I\.u2 and liu?, mono- 
tonically increase with ji. On the other hand, the peak values of the base-flow-independent 
power spectral densities H^i and lirj-, = ^, j = 2, 3}, first increase with an increase 
in /i and then plateau after ji becomes sufficiently large. The monotonic increase of Tiu2 
and n^3 with the elasticity number demonstrates the significance of the coupling opera- 
tor Cp2, which captures the work done by the polymer stresses on the flow. An in-depth 
study of the physical mechanisms behind this viscoelastic amplification is given in § |6] 
Third, the spanwise wavenumbers corresponding to the maxima in the components of the 
frequency response increase with an increase in ji and a decrease in (5. This suggests that 
the dominant structures become less spread in the spanwise direction with an increase 
in ji and a decrease in jS. 



6. Energy amplification mechanisms 

In order to elucidate the energy amplification mechanisms in Oldroyd-B fluids, we 
next analyze the Reynolds-Orr equation for streamwise-constant channel flow. As is well 
known, the Reynolds-Orr equation describes the evolution of the energy of velocity fluc- 
tuations around a given base flow condition ( Schmid fc Henningson|[2"QQl ). In our study, 
the initial conditions on velocity and polymer stress fluctuations are set to zero, but 
the flow is driven by the spatio-temporal stochastic body forcing d. This random body 
forcing generates the velocity field v and polymer stresses r, which are also of stochastic 
nature ( Farrell fc Ioannou||1993 ). 
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Figure 6. Variation in Qmax with /i for Couette (first row) and Poiseuille (second row) flows 

with /3 = {0.1, 0.5, 0.9}. 



The energy-evolution equation is derived by multiplying the Navier-Stokes equations by 
the velocity vector, followed by integration over the wall- normal direction and ensemble 
averaging in time. The equations are further simplified using the divergence theorem and 
the boundary conditions on v. For streamwise-constant perturbations, the Reynolds-Orr 
equation for an Oldroyd-B fluid is given by: 



IdE Pkl ^ , rrf X / ^\ 

2dF = -^^-^^'^^^ + ^^'^^ + ^e 
1 - p 



E - {U.U'V) + (V,d) + -^{{u.dyyU) + {v.dyyV) + {w , dyyW)) 



Re 



{{Txy.dyU) + (Tyz.dyV) + {Txz.dyW) + l/c^ ( (t^^ , ii) + {Tyz.v) + {Tzz.w))). 



where (•, •) denotes integration in y and ensemble averaging in t, and E = E{t^ kz]Re^ /i) 
represents the kinetic energy, that is 

E = (v,v) = J ^S{v''{y,k,,t)v{y,k,,t)) dy. 

Here, the asterisk denotes the complex-conjugate-transpose of vector v, and E denotes 
ensemble averaging (McComblpQlj ) 



f (v(-,t)) = hm 1 / v(-,t + r)dr. 
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In the steady- state limit, the kinetic energy E is given by 
Re 1 



/5 



{{Txy.dyu) + {Tyz.dyv) + {T^z^dyw) + l/c^ ( (t^^ , ix) + (r^^^,^) + (r^^,^))). 



In this equation, the first two terms on the right-hand-side are contributions due to the 
Reynolds stress and the work done by the body forces; the third group of terms accounts 
for viscous dissipation, and the last group of terms corresponds to the work done by the 
polymer stresses on the flow. This last contribution to kinetic energy is - in general - 
sign indefinite and it is also referred to as an energy-exchange term ( Doering et a/. [2006 ) . 
Two points are worth highlighting. First, as in Newtonian fluids (iSchmid & Henningson 



2001 ), the nonlinear terms do not contribute to the kinetic energy (they are conservative 
and only redistribute energy between different modes). Thus, the Reynolds-Orr equations 
for viscoelastic fluids derived using linearized and fully non-linear equations correspond 
to each other. Second, the steady-state energy E defined here is exactly the same as 



the ensemble- average energy density defined in Hoda et al. ( 2QQ8p . As a matter of fact, 
this quantity is precisely determined by the H2 norm of a stochastically forced linearized 
system. Thus, the energy density of the LNSE can be efficiently computed using Lya- 
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punov equations, which circumvents the need for running costly stochastic numerical 
simulations. 

The explicit scaling of the steady-state energy density with Re is given in § [4l Based 
on this, for asymptotically large times, E can be written as: 

E(kz]Re,(3,ii) = f{kz;f3,fi)Re + g{kz; f3, i^)Re^ , 



where functions / and g correspond to the Reynolds-number-independent terms in Eq. (4.3 ) 
This equation suggests that at higher Re values, g is expected to contribute most to the 
energy, whereas at smaller Re values, / is expected to contribute most to the energy; 
we will show that the latter observation holds only at moderate values of /i and p. From 
§[4j it follows that only streamwise velocity contributes to g; this contribution arises due 
to amplification from d2 and ds to u. On the other hand, amplification from di to u 
and {d2,ds) to (v^w) is captured by function /. The explicit i?e-scaling of various terms 
in the steady-state Reynolds- Orr equation can be obtained by analyzing the Lyapunov 
equation (see Appendix [B|; let the overbar designate the velocity and polymer stress 
fluctuations at i?e = 1, and let ui and U2^3 denote the components of u arising due to 
the action of di and {d2^ds)^ respectively. Then, the i?e-independent terms f{kz]/3,fi) 
and g{kz] /i) in the expression for the steady-state energy are given by 

/ = ^(^^^'^) Pii^l^^yyUl) + {V.dyyV) + {W , dyyW)) " {l - P) { {f^y , OyUl) 

+ {fyz.dyV) + {fxz.dyW) + lkz{{fxz,Ul) + {fyz.v) + (f^;, , iD ) ) ) ^ , 



3 = ^(~(^'^'^) + P {u2,3,dyyU2,3) - (1 - P){{rxy,dyU2,3) + (^a;^ , ^2,3) )) 



1 

Pk^z 

We note that function / is the same for all channel flows and only function g depends 
on the underlying base flow; for further analysis of function ^, we restrict our attention 
to Couette flow. Unless noted otherwise, all plots in this section are given in the log-log 
scale. Figures [S] and [9] show the variations in / and g with kz at different p and /i values. 
Below, we discuss the important observations concerning the results. 

We note that the magnitudes of both / and g increase with an increase in ji and a 
decrease in p. This indicates that energy amplification becomes weaker as one approaches 
the Newtonian fluid limit. Furthermore, in all the cases, / monotonically decreases with 
kz^ while g achieves a maximum at 0(1) values of kz. This suggests that the contri- 
bution fr om g is responsibl e for the energy density peaks observed at higher Reynolds 
numbers ( Hoda et fl/.||2Q Q8'), confirming our earlier claim that at higher Re, E ^ gRe^ . 

Plots in figure |9] suggest that g increases monotonically with /i, while / reaches a 
saturation limit for sufficiently large values of ji. Thus, even in low inertial regimes, 
the contribution of g to the energy density can be significant if the elasticity number, 
/i, is large enough. In particular, this demonstrates that energy density peaks observed 



by |Hoda et al\ ( |2QQ8[ ) at {k^ = 0, Re = 0.1, p = 0.1, /i = 10^} arise due to the 
contribution of g. 

Moreover, function g approximately scales linearly with /i for large values of elasticity 
number. Figure |9] shows the /C;2-dependence of g{kz', /3, in Couette flow with p = 0.1 
for five different values of the elasticity number, ja = {0.1, 1, 10, 10^, 10^}. It is evident 
that the five curves almost collapse onto each other. This is a remarkable discovery in 
view of rather complicated dependence of the underlying equations on /i and the range of 
elasticity numbers considered. Our ongoing theoretical effort is directed towards develop- 
ment of an explicit scaling of / and g with /i; we conjecture that - for large enough values 
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of elasticity number - / approximately becomes /i-independent, while g approximately 
scales linearly with ja. This would suggest the following approximate scaling of the energy 
density with elasticity number for /i 1, E{kz]Re^ ji) ^ Ref{kz] P) + iJiRe^g{kz] 
where / and g are /i-independent functions. If this scaling turns out to be correct, our 
results would indicate an interesting interplay between inertial and viscoelastic effects in 



energy amplification of Oldroyd-B fluids with low Reynolds/high elasticity numbers (Lar- 
SQii||2Q00| [Groisman fc Steinberg|[2000D . 



We have also closely examined the terms in the steady-state Reynolds-Orr equation 



that contribute to the i?e^-scaling of energy (results not shown; for details see Hoda 
(2008)). We observed that the Reynolds-stress and the energy-exchange terms, g(u,U'v) 
and g{T^y^dyu)i cire positive, suggesting that they lead to energy amplification; the vis- 
cous dissipation and the energy-exchange terms, g{u,dyyu) ^tnd g{T^^;ikzu)i ^ire negative, 
suggesting that they lead to energy suppression. Furthermore, the absolute values of all 
the terms contributing to g increase with an increase in fi and a decrease in [5. The 
Reynolds stress term, g{u,U'v)i appears to reach a saturation limit for sufficiently large 
while there does not appear to be any upper bound on the absolute values of the 
other terms contributing to i?e^-scaling as ji increases. In view of this, we conclude that 
the energy-exchange terms are mainly responsible for energy amplification in flows with 
pronounced elasticity effects. 



Figure 10 shows the /c^-dependence of the energy-exchange term g{T^y^dyu) + 9{T^zAk^u) 
in Couette now with f3 = 0.1 and /i = {0.1, 1, 10, 10^, 10^}. This term contributes to the 
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Figure 10. The /c^ -dependence of the energy-exchange term 9(^^^^ q^^) + 9{Ta,z,ikzu) in Couette 

flow with f3 = 0.1, /i = {0.1, 1, 10, 10^, 10^}. The left plot is shown in the log-linear scale and 
the right plot is shown in the log-log scale. 



i?e^-scaling of the steady-state energy and it can lead to energy suppression (for small ja) 
or energy amplification (for large ja). Also, the left plot illustrates the sign-indefiniteness 
of the energy-exchange term for moderate values of /i. Clearly, the role of the energy- 
exchange term becomes more prominent with increase in elasticity number; in particular, 
for /i 1 this term creates a much larger contribution to the energy amplification than 
the Reynolds stress term. 



It is interesting to contrast the above results with those of Sadanandan fc Sureshku- 



( 2002 ) , who studied budgets of the perturbation vorticity and kinetic energy in plane 
Poiseuille flow of Oldroyd-B fluids. Their analysis made use of the linearized equations 
and a normal- mode decomposition, and considered two-dimensional perturbations. They 
observed that the temporal frequency near the critical conditions is a non-monotonic 
function of viscosity ratio, and that perturbation shear stresses have a destabilizing 
effect. In our work, which does not use a normal- mode decomposition and considers 
three-dimensional streamwise constant perturbations, l^max is observed to be a mono- 
tonic function of P (but a non-monotonic function of /i), and perturbation shear stresses 
are found to create larger energy amplification. We ascribe the differences in these ob- 
servations to the different assumptions used in the two studies. 

Finally, we briefly comment on the terms that contribute to the i?e-scaling of energy. 
Apart from work done by the body forces and —(1 — /3)/{/3kl) {ryz, dyw)^ all other terms 
contributing to / are negative, which suggests that they suppress energy amplification; 
the most dominant of these terms are {w^dyyw)^ (ui^dyyUi)^ and (rxy^dyUi). The ab- 
solute values of all the terms, except for the work done by the body forces - which is 
independent of both P and fi - decrease with an increase in p. However, the dependence 
on /i is more subtle. The absolute values of the viscous-dissipation terms contributing 
to / increase with an increase in /i, whereas the absolute values of all other terms de- 
crease. The most striking observation is that the absolute values of both (rxy^dyUi) and 
(Txz/ikzUi) decrease with an increase in /i, whereas the absolute values of (r^;^, 1^2,3) 
and {rxz/ikzU2^3) increase with an increase in /i, suggesting competing effects of these 
energy-exchange terms at Re = 1. 
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7. Conclusions 

We have investigated the frequency responses of channel flows of Oldroyd-B fluids. 
Our analysis is based on the 2D/3C model, which simply means that only streamwise- 
constant perturbations around base values are considered. The disturbances enter the 
linearized governing equations as a body force that can vary in space and time. The 
frequency responses are described by an operator having nine blocks, with a given block 
relating the component r of the velocity perturbation to the disturbance in direction j. 
Characterization of the frequency responses can be performed through calculation of the 
Hilbert-Schmidt and H2 norms; physical interpretations of each of these were discussed. 

Our analysis shows that the frequency responses for channel flows of Oldroyd-B fluids 
scale with the Reynolds number in exactly the same way as in Newtonian fluids. Ex- 
amination of the associated block diagram reveals that the frequency responses of the 
spanwise and wall- normal disturbances to streamwise velocity scale as i?e^, whereas all 
other frequency responses scale linearly with Re. It is also seen that streamwise distur- 
bances do not affect the wall-normal and spanwise velocities. These results indicate that 
at high Reynolds numbers, channel flows of both Newtonian and viscoelastic fluids are 
most sensitive to spanwise and wall- normal disturbances, and these disturbances have 
the largest effect on the streamwise velocity. This is also reflected in the behavior of the 
various norms, which have the same scalings with Re as is seen for Newtonian fluids. The 
presence of viscoelasticity causes the frequency responses and their norms to depend on 
two additional parameters: the elasticity number and viscosity ratio. 

A parametric study of the Hilbert-Schmidt norms, or power spectral densities, shows 
that these peak at non-zero temporal frequencies for viscoelastic fluids. For Newtonian 
fluids, the peaks occur at zero frequency, indicating that viscoelasticity reduces the time 
scales over which disturbances develop. The frequencies at which the peaks occur decrease 
as the solvent contribution to the total viscosity increases, and the peaks are only present 
if the elasticity number is greater than a critical value. An analytical expression which 
was derived describes both of these trends accurately. The peak values of the power 
spectral densities and the corresponding spanwise wavenumbers increase with an increase 
in elasticity number and a decrease in viscosity ratio. This indicates that elasticity helps 
amplify spanwise-dependent disturbances and reduces their length scale. 

Since each component of the frequency response has its own peak in the Hilbert- 
Schmidt and H2 norms, different ampliflcation mechanisms may occur in a given flow. 
The componentwise analysis used in the present work makes it easy to detect this, whereas 



the aggregate analysis used in our previous work does not ( Hoda et a/. 2008) . The study 
of the power spectral and steady-state energy densities also shows that at low Reynolds 
and elasticity numbers, the flow is most sensitive to streamwise disturbances, and these 
affect the streamwise velocity the most. At low Reynolds and high elasticity numbers, the 
flow is most sensitive to wall- normal and spanwise disturbances, and these again affect 
the streamwise velocity the most. 

To elucidate how elasticity amplifles disturbances, we have analyzed a Reynolds-Orr 
equation for the kinetic energy. Its steady-state value is exactly the same as the H2 
norm. Our analysis shows that the energy-exchange term involving the polymer stress 
component Txy and the wall- normal gradient of the streamwise velocity dyU promotes 
energy ampliflcation. In contrast, the energy-exchange term involving the polymer stress 
component and the spanwise gradient of the streamwise velocity ikzU leads to en- 
ergy suppression. The energy-exchange term that promotes energy ampliflcation becomes 
increasingly important relative to the Reynolds stress term as the elasticity number in- 
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creases, and is thus the main driving force for amphfication in flows with strong vis- 
coelastic effects. 

The results of the present work significantly extend our knowledge of how viscoelastic 
channel flows respond to external disturbances. They clearly demonstrate how streamwise- 
constant disturbances can produce significant energy amplification, and clarify the re- 
lationship between the components of the disturbances and velocity field. Notably, the 
results cover both inertia- and elasticity-dominated flows. We expect that they will pro- 
vide a useful basis for studies exploring nonlinear aspects of transition to inertial and/or 
elastic turbulence in channel flows of viscoelastic fluids. It might also be possible to extend 
the ideas developed here to study interfacial instabilities, where the linearized governing 
equations are also non- normal in nature (Lin, Khomami fc Sureshkumar^^2QQ4) . 
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Appendix A. Frequency response operators 

The Reynolds-number-independent frequency response operators H^j can be obtained 
by applying the temporal Fourier transform to Eqs. (2.3). Equation (2.3b) can be used to 

= V 

(Al) 



express polymer stresses 02 = [ryy Tyz Tzz]^ in terms of wall- normal velocity (j)! = v 

F21 



02 = 



1 + [jj^n 



Q — uoRe. 



Substitution of Eq. (Al) into the Fourier transform of Eq. (2.3a) yields 

-1 



1 



F12F21 (B2(i2 + Bs^s) . 



(A 2) 



Using Eq. (2.3/) and the fact that F12F21 = Fn = A~ A , the expressions for opera- 
tors H^j(/c2, 1]; /i), {r = v^w] j = 2,3}, are readily established (cf. §[4|. 

The following relation between the wall- normal velocity /vorticity {di^J^s) and polymer 
stresses 04 = [ Txy Txz ] can be established by substituting Eq. ( A 1 ) to the Fourier 
transform of Eq. (2.3d) 
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Substituting the above equation in the Fourier transform of Eq. ( |2.3c ) yields 
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Now, since u = C^^a, the expressions for operators fLujikz, /5, /i), {j = 1, 2, 3}, given 
in § |4] follow immediately from this equation and the fact that in Couette and Poiseuille 
flows we have {F34F43 = F33 = A, F34F41 = 0, F34F42F21 = ik, {U\y)A + 2U'\y)dy)}. 
The latter expression appears in the coupling operator Gp2 of § |4] A careful examina- 
tion of the governing equations shows that this term arises due to the tensor involving 
polymer stress fluctuations and gradients in the base velocity proflle (the r • Vv term in 
Eq. (I2II)). 



Appendix B. Determination of the H2 norm using Lyapunov equation 

Next, we outline the procedure that is most convenient for determination of the 
Reynolds-number-independent functions / and g in the expression for the energy am- 
pliflcation. This method exploits the fact that the steady- state energy density of the 
frequency response operator H can be determined from the Lyapunov equation, and 
it is very suitable for uncovering the explicit scaling with Re of various terms in the 
steady-state Reynolds-Orr equation. 
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B.l. Lyapunov equation and steady-state energy density 
Equations (2.3a)-(2.3d) from §[2]can be compactly rewritten as 
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Furthermore, the structure of operators Bn and B21 yields 
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where the asterisk denotes the adjoint of a given operator. From the deflnitions of oper- 
ators Bj, j = 1,2,3, and their respective adjoints ( Jovanovic fc Bamieh||2005 ) we have 
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The partition of operators Mn and M22 is done conformably with the partition of vectors 



^1 



^1 02^ 



1 T 



and = 



, respectively. 

The steady-state energy density (that is, the H2 norm of operator H) can be expressed 
in terms of the solution to the following operator Lyapunov equation 



AP + PA* 



-M 



as E{kz] Re, (3, jii) = trace (NP), where P and N are 2x2 self-adjoint block operators 
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with the fohowing structure 
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Using the structure of operators P and N, it fohows that E{kz] Re^ /3, /i) is determined by 
E{kz]Re^ 13^ 11) = trace (NP) = trace (NuPn) + trace (N22P22)- The same argument 
as in |Bamieh fc Dahieh| ( |2QQl| ); |Jovanovic fc Bamieh (|2QQ5) can be used to show {Pn = 
Re X, P21 = i?e^ Y, P22 = i?e Zi + i?e^ Z2}, where the Reynolds- number- independent 
operators X, Y, and Z satisfy the fohowing system of conveniently coupled equations 
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Hence, 
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= (trace (NiiX) 
= Ref{kz',P,fi) + 



+ trace (N22Z1)) 
Re^ g(kz]p,ii), 



Re^ trace (N22Z2) 



which provides an efficient way of computing functions / and g in the expression for the 
steady-state energy density. 

B.2. Lyapunov equation and steady-state Reynolds-Orr equation 

We next illustrate how the terms on the right-hand side of the steady-state Reynolds- 
Orr equation can be obtained from the solution of the Lyapunov equation. This allows 
us to determine the explicit i?e-scaling of these terms, which clarifies the importance of 
different energy amplification mechanisms. 

We begin by observing that the solution of Lyapunov equation, P, represents the 
steady-state correlation operator of = [ Lp{ j 5 ^^^^ 



L ' xy ' xz J 

kx — 0, i/- = C-j 

^7 "^xy ■) cind Txz 



P = lim £{>p{;t)<E>>p{;t)}, 

t ^ 00 

where cp ^ cp denotes the tensor product of (p with itself. This operator contains all 
the second-order steady-state statistics of the velocity and polymer stress fields. For 
example, operator P22 carries information about the steady-state correlation of cp = 
with itself, and a simple kinematic relationship between cOy and u (at 
ujy = —{i/kz)u;y) can be used to obtain all possible correlations between 
In particular, owing to the convenient scaling of P22 with the Reynolds 
number (P22 = ReZi + Re^ Z2), we conclude that all these correlations contain a part 
that scales as Re and a part that scales as Re^. (The streamwise force is responsible for the 
i?e-scaling and the wall- normal and spanwise forces are responsible for the i?e^-scaling.) 
Now, since the inner product of two fields is equal to the trace of their outer (tensor) 
product, we conclude that the terms (u^dyyu)^ {Txy^dyu)^ and (rxz^^kzu) in the steady- 
state Reynolds-Orr equation scale as fi{kz; fi)Re + gi{kz; /J^)Re^ ^ {i = (u^dyyu), 
i^xy^dyu), {Txz/ikzu)}. Since the Reynolds stress term —{u^U'v) can be determined 
from P21 = i^e^Y, it follows that its contribution to the steady-state energy scales as 
Re^ (cf. §[6|. The contributions of all other terms scale as Re. 



